Canonical Expansion of VT— Symmetric Operators 

and Perturbation Theory 



E.Caliceti* 

Dipartimento di Matematica, Universita di Bologna 
40127 Bologna, Italy 
S.Graffit 

Department of Mathematics and Computer Science 
Emory University, Atlanta, Ga 30322. U.S.A. 



Abstract 

Let H be any VT symmetric Schrodinger operator of the type — h 2 A + (xi + . . .+ 
x d) + igW(xi, . . . , Xd) on L 2 (R rf ), where W is any odd homogeneous polynomial 
and §eR. It is proved that VH is self-adjoint and that its eigenvalues coincide 
(up to a sign) with the singular values of H, i.e. the eigenvalues of V H*H. 
Moreover we explicitly construct the canonical expansion of H and determine 
the singular values [ij of H through the Borel summability of their divergent 
perturbation theory. The singular values yield estimates of the location of the 
eigenvalues Xj of H by Weyl's inequalities. 



1 Introduction and statement of the results 

A Schrodinger operator H = — A + V acting onH = L 2 (R d ) is called VT— symmetric 
if it is left invariant by the VT operation. While generally speaking V could be the 
parity operator with respect to at least one variable, here for the sake of simplicity we 
consider only the case in which V is the parity operator with respect to all variables, 
(Vu)(xi, . . . , Xd) = u(—xi,...,—Xd), and T the complex conjugation (equivalent to 
time-reversal symmetry) (Tu)(xi, . . . , Xd) '■= u(xx, . . . , Xd)- The condition 

V(-x u . . . , -x d ) = V{x u ...,x d ) (1.1) 
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defines the VT— symmetry on the potential V(x±, . . . ,Xd). The PT-symmetric oper- 
ators are currently the object of intense investigation because, while not self-adjoint, 
they admit in many circumstances a real spectrum. Hence the investigation is moti- 
vated (at least partially), by an attempt to remove the self-adjointess condition on the 
observables of standard quantum mechanics (see e.g. [1], [2], [3], [4], [5], [6], [7], [8], [9]). 

The simplest and most studied class of VT symmetric operators is represented 
by the odd anharmonic oscillators with purely imaginary coupling in dimension one, 
namely the maximal differential operators in L 2 (R) 

d 2 

Hu(x) :=[-— + x 2 + igx 2m+1 l g G R, m = l,2,... (1.2) 
dx l 

It has long been conjectured (Bessis Zinn- Justin), and recently proved [10], [11], that 
the spectrum a(H) is real for all g; there are however examples of one-dimensional 
"PT-symmetric operators with complex eigenvalues [5]. 

Now recall that there is a natural additional notion of spectrum associated with a 
non-normal operator T in a Hilbert space which is by construction real. Any closed 
operator T admits a polar decomposition ([12], Chapt. VI. 7) T = U\T\, where \T\ is 
self-adjoint and U is unitary. The modulus of T is the self-adjoint operator \T\ = VT*T. 
The (obviously real and positive) eigenvalues of |T| are called the singular values of 
T. In this paper we consider the self-adjoint operator \J H*H; its eigenvalues j = 
0, 1, . . ., necessarily real and positive, are the by definition the singular values of H . A 
first immediate question arising in this context is to determine how these singular values 
are related to the PT-symmetry of H. A related question is the explicit construction of 
the canonical expansion of H (see e.g. [12]) in terms of the spectral analysis of \/H*H, 
which entails the diagonalization of H with respect to a pair of dual bases (which do 
not form a biorthogonal pair); a further one is the actual computation of the singular 
values. The determination of the singular values reflects directly on the object of 
physical interest, namely the eigenvalues A^; j = 0, 1, ... of H . If the eigenvalues and 
the singular values are ordered according to increasing modulus, the Weyl inequalities 
(see e.g. [18]) indeed yield 
k k 

X!l A jl^X^i> | Ai - - - Afc| < a*i - - - a**., A; = 1,2,... (1.3) 

i=i i=i 
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We intend in this paper to give a reply to these questions for the most general class 
of odd anharmonic oscillator in R d . Namely, we consider in L 2 (R, d ) the Schrodinger 
operator family 

H(g)u(x) := H u(x) + igW(x)u(x), x = (x u . . . , x d ) G R d (1.4) 

Here: 

1. W is a real homogenous polynomial of odd order 2K +1, K = 1,2 . . .; 

W(Xx) = X 2K+1 W(x) 

2. H is the Schrodinger operator of the harmonic oscillator in R d : 

H u(x) = —Au(x) + x 2 u(x), x 2 := x\ + . . . + x 2 d (1.5) 

Under these conditions the operator family H(g), which is obviously PT-symmetric 
(see below for the mathematical definition), but non self-adjoint, enjoys the following 
properties (proved in [13] for d — 1 and in [14] for d > 1; see below for a more detailed 
statement): 

1. The operator H(g), defined as the closure of the minimal differential operator 
H(g)u = —Au(x) + x 2 u(x) + igW(x), u G C£°(R d ), generates a holomorphic 
operator family with compact resolvents with respect to g in some domain S C C, 
with H(g)* = H(g). An operator family T(g) depending on the complex variable 
g G f2, where f2 C C is open is holomorphic (see [12], VII. 1) if the scalar products 
(u,T(g)v) are holomorphic functions of g G V (u, v) G T(g) and the resolvent 
[T(g) — zl]~ l exist for at least one g G fl 

2. All eigenvalues of H := H(0) are stable with respect to the operator family H(g). 
This means (see e.g. [12], VIII. 1) that if Ao is any eigenvalue of H(0) of multiplicity 
m, there is B(X ) > such that H(g) has exactly m (repeated) eigenvalues 
^'(fiO : J ' = 1? • • • ? m near -^o f° r 9 ^ <S, \g\ < B(Xq), and lim Xj(g) = X . 

3. The (Rayleigh-Schrodinger) perturbation series of each eigenvalue X(g) of H{g) 
is Borel summable to X(g). 
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We denote (J,j(g) : j = 0,1,... the singular values of H(g) : g G R, i.e. the 
eigenvalues of ^H(g)*H(g) = ^H(-g)H(g). 

Our first result concerns the identification of the singular values as the eigenvalues of 
a self-adjoint operator directly associated with H{g) by the operator-theoretic imple- 
mentation of the recently isolated pseudohermiticity notion ([1], [15], [16], [17]) in terms 
of the V symmetry itself. 

Consider indeed the operator family Q(g) :— VH(g). We will show that D(Q(g)) = 
D(H(g)). Since [V, H ] = 0, the explicit action of Q(g) is 

Q(g)u(x) = H u(—x)+igW(—x)u( y —x) 

= H u(-x) - igW{x)u{-x) = H(-g)Vu(x) 

Then we have: 

Theorem 1.1 Let Q(g) be defined as above and Q'(g) := H(g)V. Then: 

1. If g G R the operator families Q(g) and Q'(g) are self-adjoint. 

2. The operator family Q(g) defined on D(H(g)) is holomorphic with compact re- 
solvents at least for g in a neighbourhood of R + . 

3. If g G R the eigenvalues of Q(g) and of ^ H(g)*H(g) coincide (up to the sign); 
Remarks 

1. H(g)* = H(—g) for g G R by "PT-symmetry. Hence the relation Q(g) = 
VH(g) = H(-g)V = H(g)*V can be equivalently written VE{g)V' 1 = H(g)* 
which is the V— pseudohermiticity property of H(g) [16]. 

2. The eigenvalues \x of the operator Q(g) clearly solve the generalized spectral 
problem H(g)u = \xVu (for this notion, see [12], SVII.6). Explicitly: 

(H + igW)u{x) = fi{Vu){x) (1.6) 

By the above theorem the singular values coincide (up to a sign) with the gener- 
alized eigenvalues. 
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As a consequence of this, we obtain the explicit canonical expansion of H(g) in terms 
of the eigenvectors ip k of Q and of the V operation: 

Corollary 1.2 Let {ip k (g)} '■ k — 0,1, ... be the eigenvectors ofQ(g), and [i k the corre- 
sponding eigenvalues (counting multiplicy). Then H(g) admits the following canonical 
expansion 

oo 

H{g)u = £ fi k (u, 4> k )P4> k , u E D(H(g)) (1.7) 

fc=0 

Remarks 

1. Since Vi/j n (x) = ip n (—x), the canonical expansion (1.7) entails that H(g) can be 
diagonalized in terms of the (repeated) real singular values //„ and of the pair of 
orthonormal bases {ip n (x)} and {ip n (— x)}. 

2. For a general operator with compact resolvent the canonical expansion reads 

oo 

Tu = y £»k{u,'<l>M, ueD(T) (1.8) 

fc=0 

Here {fi k } is the sequence of singular values of T, i\) k the corresponding eigenvec- 
tors, but the dual basis {ip' k } is a priori unknown. In this case it is simply the 
P-dual basis Vip k . Remark that the orthogonal bases ip k and Vip k do not form a 
biorthogonal set. 

3. The expansion (1.7) is useful even when all eigenvalues of H(g) are real, because 
H(g) is not normal and the spectral theorem does not hold. 

4. Finally we note the following relation involving nonzero eigenvalues and eigen- 
vectors on one side and nonzero singular values and corresponding eigenvectors 
on the other side: if H(g)ip k = /i k (g)Vip k , and H(g)<f>i = Xi(g)(f>i, then: 

\i(g){<p h VMg)) = ^(g){<PuMg)) (i-9) 



One has indeed (omitting the (/-dependence): 
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Our third result deals with the actual computation of the singular values Hj{g). To 
formulate the result, remark that the closed subspaces VTC and (1 —V)H are invariant 
under H because [V, H ] = . The operator VH has the same eigenvectors of H Q , 
but the eigenvalues A; = 2/ x + . . . + 2l d + d, l k = 0, 1, . . ., k = 1, . . . , d, of H split 
into even and odd eigenvalues. More precisely, introduce the usual principal quantum 

number I :— l\ + . . . + : I — 0, 1, Then the eigenvalues of H are A/ = 21 + d, with 

multiplicity m(l) = l d ~ x . The eigenvalues of VHq are 

[ -(21 + d), I odd v ' 

The corresponding eigenvectors will be V even and V odd, respectively. We then have: 

Theorem 1.3 

1. All eigenvalues A/ of VH are stable as eigenvalues fj,j(g) : j = l,...,m(l) of 
Q(g) as \g\ — > 0, <? G <Si U <5>2 where: 

:= {^GC\{0}:-^<arg^<|} (1.11) 

5 2 := {^GC\{0}:|<arg^<^} (1.12) 

2. All eigenvalues (J>j(g) : j — 1, . . . ,m(l) are holomorphic on the Riemann surface 
sector 

S K ,s :={geC:0<\g\< B(5); -{2K + 1)^ + 5 < arg (<?) < (2tf + 1)^ - 5} 
where 5 > is arbitrary. 

3. The Rayleigh-Schrddinger perturbation expansion for any eigenvalue ^j(g) : j = 
1, . . . , m(Z) ofQ(g) near the eigenvalue A; ofVH for \g\ small is Borel summable 
to Hj{g) :j = l,...,m(l). 

Remark 

Let n(g) be a singular value near an unperturbed eigenvalue A. The Borel summability 

(see e.g. [19], Chapter XII. 5) means that it can be uniquely reconstructed through its 

00 

divergent perturbation expansion fJ> s g s , f>o = A in the following way: 

s=0 

1 r°° 1 / / 

M = - Mg^e- 1 t-^dt (1.13) 

C£ J 
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Here q = and {13(9), the Borel transform of order q of the perturbation series, 

is defined by the power series 

M9)= 5frjtt)j 9 * 

which has a positive radius of convergence. The proof of (1.13) consists precisely in 
showing that fiB(g) has analytic continuation along the real positive axis and that the 
integral converges for some < g < B, B > 0. 
Example 

The Henon-Heiles potential, i.e. the third degree polynomial in R 2 

W(x) = x\x 2 

2 Proof of the results 

Let us begin by a more detailed quotation of Theorem 1.1 of [14]. The results are more 
conveniently formulated in the variable j3 — ig instead of g. 

Let (3 G C, < |arg/3| < 7r, and let H{f3) denote the minimal differential operator 
in L 2 {R d ) defined by -A + x 2 + f3W{x) on C °°(R d ), with x 2 = xf + . . . + x\. Then 

(Nl) H(f3) is closable. Denote H(f3) its closure. 

(N2) H(f3) represents a pair of type-A holomorphic families in the sense of Kato for 
< arg/? < 7r and — ir < arg/3 < 0, respectively, with H(f3)* = H(f3). Recall that 
an operator family T(g) depending on the complex variable g belonging to some 
open set Q C C is called type-A holomorphic if its domain D does not depend 
on g and the scalar products (u,T(g)) are holomorphic functions for g G D 
V (u,v) G D. A general theorem of Kato ([12], VII. 2) states that the isolated 
eigenvalues of a type-A holomorphic family are locally holomorphic functions of 
g G D with at most algebraic branch points. 

(N3) H{(3) has compact resolvent V/3 G C, < |arg/3| < it. 

(N4) All eigenvalues of H = H(0) are stable with respect to the operator family H(j3) 
for (3 — > 0, < |arg/3| < 7r; 
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(N5) Let P E C, a E C, < |arg/?| < n, -it + arg/3 < arga < arg/3, and let H a {(3) 
denote the minimal differential operator in L 2 (R d ) defined by — A + ax 2 + /3W(x) 
on C^(~R d ). Then H a (f3) is sectorial (and hence closable) because its numerical 
range is contained in the half-plane {z E C : — ir + arg/3 < arga < arg/3}; 

(N6) Let H a (P) denote the closure of H a ((5). Let a E C, a - oo,0]. Then the 
operator family j3 \— > H a ((3) is type- A holomorphic with compact resolvents for 
/3 G C CT := {/3 E C : < arg/3 - argcr < vr}. Moreover if (3 E C,Im/3 > 0, the 
operator family a \— > H a {(3) is type- A holomorphic with compact resolvents in 
the half-plane = {cr G C : < arg/3 — argc < it} 

Let us now introduce the operator 

H((3, 6) = e- 2d A + e 29 x 2 + (3e 2 ^ d W{x) := e^K^, 6) (2.1) 

For 9 E R H(P,9) is unitarily equivalent to H(/3), Im/5 > 0, via the dilation operator 
defined by 

(U(9)4>)(x) = e de ^(e e x), Wip E L 2 (R d ) (2.2) 

As a consequence of (N6) (see again [14], or also [20], where all details are worked out 
for d — 1, and where the reader is referred also for the proof of statement ( N8) below) 
we have: 

(N7) H(j3,9) defined on D(H(j3)) represents a type-A holomorphic family with com- 
pact resolvents for (5 and 9 such that s = arg/5, t = lm9 are variable in the 
parallelogram 1Z defined as 

K = {(s, t) E R 2 : < (2K - l)t + s < ir, < (2K + 3)t + s < vr}, (2.3) 

Moreover is a core of H(/3, 9). The spectrum of H(j3, 9) does not depend on 
9. Note that (s,t) E 1Z entails that the maximal range of (5 is —{2K — 1)tt/A < 
arg/3 < (2K + 3)rr/4 and that the maximal range of 9 is — tt/A < lm9 < ir/4; 

N8) Let and 9 be such that (s, t) G K. Then: 
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(i) If A i a(K(0, 9)), then A G A, where: 

K:={zeC:zi a(K(f3,9)); \\ [z - (K(P,9)]- l \\ (2.4) 
is uniformly bounded for ^ 0}; 

(ii) If A G a(K(0, 9)), then A is stable with respect to the operator family 
K((3,9). 

(N7) and (N8) entail: 

(N9) Let (3 G C with < arg(/3) < n. Then for any 5 > and any eigenvalue X(g) of 
H(f3) there exists p > such that the function A(/3), a priori holomorphic for < 
\g\ < p, 5 < arg(/3) < tt — 5, has an analytic continuation to the Riemann surface 
sector S KjS := {(3 G C : < \(3\ < p; -(2K - 1)^ + 5 < arg(/3) < (2K + 3)^-5}. 

Remarks 

1. The stability statement means the following: if r > is sufficiently small, so that 
the only eigenvalue of K(0, 9) enclosed in r r := {z G C : \z — A| = r} is A, then 
there is B > such that for \/3\ < B dimP(/3, 9) = dimP(0, 9), where 

P(P, 9 ) = ^-j^ z - ^ ( 2 - 5 ) 

is the spectral projection of K(f3, 9) corresponding to the points of the spectrum 
enclosed in Y r C C \ a(K((3, 9)). Similarly for P(0, 9). 

2. Starting from the operator H(f3), Im/5 < 0, analogous results hold for the operator 
family H([3, 9) where this time (3 and 9 are such that s = arg/3 and t = lm9 
describe the parallelogram 

K 1 = {(s, t) G R 2 : -7T < (2K - l)t + s < 0, -7T < (2K + 3)t + s < 0}. (2.6) 

Moreover, H((3,9)* = H(p,9). 

We now set (3 = ig and with slight abuse of notation the operator H((3) = H{ig) will 
be denoted H(g). 
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Let once again V denote the parity operator in H 

Vi/j(x) = ij)(-x), G H 

P is a self-adjoint, unitary involution, i.e. V 2 = I, and VW(x) = —W(x)\/x G H d . 
To show Theorem 1.1, let us first state and prove the following preliminary result: 

Proposition 2.1 Let Si, S 2 be the complex sectors defined by (1.11,1.12). Then: 

(1) D{VH{g)) = D{H{g)V) = D{H{g)*V) = D(VH(g)*) = D(H(g)) for all g G 
Si U S 2 ; 

(2) VH(g)) = H{-g)V for all g E S 1 L)S 2 . In particular, for g E R, VH{g) = 
H(g)*V whence VH(g)V = H{g)*, i.e. H{g) and H(g)* are unitarily equivalent; 

(3) VH{g)i> = H(g)P$, G D{H{g)), Wg G R. 
Proof 

(1) Since H(g)* = H(—g), and D(H(g)) is independent of g G Si U «S 2 , it is enough to 
prove that, for all g G S\ U S 2 : 

(a) D(VH(g)) = D(H(g)); (b) D(H(g)V) = D(H(g)) V g G S 1 U S 2 

(a) follows from D(V) = H. As for (b) notice that u G D(H(g)) if and only if 
3{m„} G C^°(R d ) such that u n — > u and H(g)u n — > t> = H(g)u. Then w„(— x) G 
C °°(R d ) -> u(-a:) and H(-g)u n (-x) -> u(-a:). Thus = u(-ar) G D(H(-g)) = 
D(H(g)), i.e. « G D(H(g)V). Conversely, if u G D(H(g)V) then u(-a;) G D(H(g)) 
and u G D(H(-g)) = D(H(g)), whence D{H(g)) = D(H(g)V). 

(2) From (1) we have D(VH(g)) = D(H(-g)V) = D(H(g)); moreover C^R^) is a 
core for both operators VH(g) and H(—g)V. Therefore it is enough to prove that 
VH(g)u = H{-g)Vu V«G C °°(R d ). Indeed, if u G C °°(R d ) then Vu G C °°(R d ) and 

VH(g)u = V(-Au + x 2 tp + igWu) = -AVu + x 2 Vu - igWVu = H{-g)Vu 

(3) Again it is enough to prove the identity for u G C^°(R d ). By direct inspection: 



VH{g)iP = H{-g)Vip = -AVip + x 2 Vip - igWVip = H{g)Vip 
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because Vip = Vip. This proves the Proposition. 
Proof of Theorem 1.1. 

1. Since V is continuous in H we have ([12], Problem 5.26) Q(g)* = (PH(g))* = 
H(g)*V = VH(g), where the last equality follows from Assertion (2) of Proposition 
2.1. Since VH(g) = Q(g), Q{g)* = Q{g)- Same argument for Q'(g). 

2. The domain of H(g) does not depend on g by (N2) for g G S\ U 1S2. Hence also 
the domain of Q(g) is g— independent. Moreover the scalar products (Q(g)u,u) are 
obviously entire holomorphic functions of g Vw G D(Q(g)). Thus Q(g) is by definition 
a type- A holomorphic family in the sense of [12] (Section VII. 1.3) in the stated domain. 
We now verify that p(Q(g)) 7^ for g belonging to a neighbourhood of R + . Since 
^ a(H ), by (2.4) with 8 = there is B > such H(g)~~ 1 is uniformly bounded in 
S := {g G 1S1U1S2, < B}. Hence \x = is not an eigenvalue of Q(g) = VH(g) because 
V is invertible. Therefore Q(g) is invertible for g G S. Now Ran(Q(g) = L 2 : if indeed 
v G L 2 , then Vv G R(H(g)) = L 2 , i.e. there exists u G D(H(g)) such that H(g)u = Vv. 
Hence VH(g))u = v and v G Ran(VH(g)). The inverse Q(g)^ 1 = H(g)~ l V is compact 
as the product of the compact operator H^g) -1 times the continuous operator V . Since 
Q{g) is self-adjoint for g G R, the compactness of the resolvent [Q(g) — z]^ 1 extends 
to all g in a neighbourhod of the real axis (see [12], Thm VII. 2. 8). 

3. Let us first prove the coincidence between the eigenvalues of Q(g) and those of 
Q\g). We have: 

VH(g)ij) = \ip <^=> H{g)ij) = XP^J {H{g)V)V^ = Q'(g)Vi/> = Wip 

Hence A is eigenvalue of Q(g) with eigenvector ip if and only if A is eigenvalue of Q'(g) 
with eigenvector Vip. 

Let now \x be any eigenvalue of Q = VH, and let ip be any corresponding eigenvector. 
Then, by the self-adjointness of VH: 

QiP = VHiP = flip =>■ H*HiP = H*VVHiP = Q 2 ^ = /i 2 ^. 

Thus fi 2 is an eigenvalue of H*H with the same eigenvector of Q. On the other hand, 
since, as we have seen, Q' 1 exists and is compact the eigenvectors of Q form a complete 
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set. Therefore fj 2 is an eigenvalue of H*H if and only if /i or —ji is an eigenvalue of Q. 
This completes the proof of Theorem 1.1. 

Proof of Corollary 1.2. 

By the spectral theorem we have, if u G D(Q): 

oo 

QU = VHU = Hn{u, ^n)^n 
n=0 

(counting multiplicities). Since VQ = H, and V is continuous: 

oo 

Hu = J2 Vniu, i>^Vi> n , V« G D(H(g)) 

n=0 

NOW (7>Vn)(z) =^„(-x) 

Define now 9) := VH((3, 9) and let us prove that this operator family enjoys the 
same properties of H(fi, 6). We have: 

Proposition 2.2 Q(f3,9) defined on D(Q(j3)) = D(H(f3)) is a type-A holomorpic fam- 
ily with compact resolvents in a neighbourhood of R + for f3 and 6 such that (s, t) G 71, 
s = arg/5 7 t = lm6. Moreover C^°(R d ) is a core of Q(f3,9). Analogous results hold for 
the operator family Q((3, 9) for f3 and 9 such that (s, t) G 7V, and Q({3, 9)* = Q((3, 9). 

Proof. 

The fact that Q(/3,0) is closed on D(H((3)) = D(H((3,9)) can be proved by the same 
argument of Proposition 2.1, (1). To complete the proof we then proceed as in Theorem 
1.1, (2). This proves the proposition. 

Proof of Theorem 1.3 

Set T(f3,0) := e 29 Q({3,9) = VK({3,9). Given the analyticity property of the operator 
family Q((3, 9), we have only to verify the analogous of N5); namely that, for all 9) 
such that (s, t) G 71, the following two properties hold: 

(i) If A £ a(T(0, 9)), then A G A 1 where: 

Ax := {z G C : z£ a(T((3,9)); \\[z - (T((3, 9)]- 1 \\ (2.7) 
is uniformly bounded for \/3 \ — > 0}; 
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(ii) If A G cr(T(0, 9)), then A is stable with respect to the operator family T(f3, 9). 

To prove these assertions, we generalize the argument of [20] valid for d = 1. First set 
p := \(3\, K(p) := K((3, 9), T(p) := T((3,9). The proof of N10) relies on the following 
results (see [21]): 

(a) \imK(p)u = K(0)u, \imK( P y U = K(0)*u, \/u G C °°(R d ); 

(b) A, ? 0; 

(c) Let x e C^°(R d ) be such that < x(%) < 1, x{%) = 1 for|a:j < 1, x{%) = 
for | a; | > 2. For h G N let Xh( x ) '■= x( x /h), and M h (x) = 1 — Then: 

(1) If p m J, and u rn G D{K{p m )) are two sequences such that ||u m || — > 1, 
M m — > weakly, and || (K(p m ))u m \\ is bounded in m, then there exists a > 
such that 

liminf ||M/j-u m || > a > 0, V/i 

m— >oo 

(2) For some z G Ai: 

lim || [Af fc , .Kr(p)] - AT(p)]- 1 1| =0 

(3) lim d h (\p) = +oo VA G C, where: 

h^>oo 
PlO 

d h (X,p) := inf{||[A - K{p)]M h u\\ : « G D(K(p)), \\M h u\\ = 1} 

Hence we must verify the analogous properties, denoted (a'), (&'), (c'), for the operator 
family T(p). Remark that, as in [20], the verification of (b') requires an argument 
completely independent of [21] because the operator family T(p) is not sectorial. We 
have: 

(a') From (a) and the continuity of V we can write 

lim T(p)u = T(0)u, \imT(pYu = T(0yu, VMC co (R d ) 

(b') First remark that G A by N9) (i) since ^ a(K(0, 9)). Then there is B > such 
that 

sup WK^ey^i < +oo. 

0<|/3|<B 
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To prove the analogous bound with T(/3, 8) in place of K(/3, 8), note that VK((3, 9)ip = 
if and only if ip — 0. Hence there exists B > such that /i — is not an eigenvalue 
of T(/3,0) for < 5. Thus T(/3, 8) is invertible. Its range is H: if u G W, then 
G Ran(K((3,8)) = H, i.e. there exists u G D(K(f3,8)) such that K((3,8)u = Vv. 
Thus VK((3,8)u = v and t> G Ran(VK((3, 8)). Finally T(/3, p) _1 = = 
Kifi^^V is uniformly bounded for < _B because P is bounded and K(f3,8)~ l is 
uniformly bounded. 

(c') Let x £ C^°(R d ) be as in (c) with the additional condition x( x ) = x{~ x )i i- e - 
V X = X- Then V Xh = Xh, and VM h = M h . We have: 

(1') Let p m I and it m G D(T(p m )) be such that ||u m || — > 1, u rn — > weakly and 
||^(Pm)«m|| < const. Vm. Then ||if(/0 m )u m || = ||T(p m )u m || < const.; hence by 
(cl) there exists a > such that 

liminf II M h u m II > a > 0, W h 

m—>oo 

(2') As proved in [21], if (c2) holds for some z G A^ then it holds for all z G A^. Thus 
we can take z = G A n Ai and we have: 

lim ||[M,,T(p)](^(p))- 1 || = lim || (M h VK(p) - VK(p)M h )(VK(p)y 1 1| 

n—*oo h— >oo 

= lim ||P[M /l ,X(p)] J ft:(p))- 1 P||=0 

ft— >oo 

where the last equality follows from the unitarity of V and (c2). 
(3') Let A G C and 

d' h (X,p) := inf{||(A - T{p))M h u\\ : u G £>(T(p)), ||M h «|| = 1} 

Then: 

||[A-T(p)]M h «|| = \\[X(l-V)+V(X-K(p))]M h u\\ > 
|| [A - tf(p)]M fc u|| - |A|||(1 - P)M^|| > || [A - K(p)]M h u\\ - |A| 

Hence d' h (X, p) > d h (X, p) — |A| and by (3) lim d' h (X, p) = +oo. The assertion 

h— >oo 

is now a direct application of [21], Theorem 5.4. This concludes the proof of 
Assertions 1 and 2 of Theorem 1.3. 
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Let us now turn to the proof of Assertion 3, i.e. the Borel summability of the eigenvalues 
of the operator family Q(g,9) := Q(f3 ) 9) for (3 = ig, —ir/A < argg < 7r/4, \g\ suitably 
small (depending on the unperturbed eigenvalue). 

To this end, we adapt to the present situation the proof [14] valid for the operator 
family H(g, 9) := H(f3 ) 9), (3 = ig, in turn based on the general argument of [22]. 

First remark that if (/3, 9) generates the parallelogram 1Z defined in (2.3) then (g, 9) 
generates the parallelogram 

K = {(s, t) E R 2 : -tt/2 < (2K -l)t + s< tt/2, -tt/2 < (2K + 3)t + s < vr/2}, (2.8) 

where now s = argg = arg/3 — tt/2. From now on, with abuse of notation, we write 
(g, 9) e TZ whenever (s, t) £ 1Z. 

Let A be an eigenvalue of H (9) := H(0,9) of multiplicity m(A) := m. Denote 
P(0, 9) the corresponding projection. By the above stability result, this means that if 
T is a circumference of radius e centered at A there is C > independent of (g,9) G 1Z 
such that, denoting Rg(z,g,9) := [Q(g,9) — z]^ 1 the resolvent of Q(g,9): 

sn P \\[Q(g,9)-z]- 1 \\ <C, \g\ - 

zero 

and that dimP(g,9) = dimP as \g\ — > 0, (p,^) G 7?., argg fixed. This time: 

P(g, e )--=^-j T R Q^ 9, 0) dz, P = P(0, 0) := ^- jT P Q (z, 0, 0) rfz (2.9) 

are the projections on the parts of a(Q(g,9)), a(VH(0,9)) enclosed in T. We recall 
that a(Q(g,9)) is independent of 9 for all ((?, 9) in the stated analyticity region, and 
that P(0,9) = P(0,9). It follows that Q(g,9) has exactly m eigenvalues (counting 
multiplicities) in T, denoted once again /ii(g) , . . . , /i m (g) . We explicitly note that, 
unlike the m — 1 case, when the unperturbed eigenvalue is degenerate, the analyticity 
of the operator family does not a priori entail the same property of the eigenvalues 
Hi(g), ■ ■ ■ ,pL m (g) : so that the analysis of [14], [22] is necessary. Following [[22], Sect. 5] 
set: 

M(g, 9) := Ran(P Q (g, 9)); D(g, 9) := P(0, 9)P(g, 9)P(0, 9) 

Under the present conditions D(g,9) is invertible on M(0) := Ran(P(0,9)). Hence 
the present problem can be reduced to a finite-dimensional one in .M(0, 9) by setting 



16 



N(g,9) := P(0,9)P(g,9)[Q(g,9)-X]P(g,9)P(0,9) 

As in [[22], Thms 4.1, 4.2] the Rayleigh-Schrodinger series for each eigenvalue /J, s (g) : 
s = 1, . . . , m near A is Borel summable upon verification of the two following assertions: 
there exist r)(5) > and a sequence of linear operators {Ei(0, 9)} in At(0, 6 1 ) such that 

(i) E(g,9) is an operator- valued analytic function for (g, 9) G 7Z; As we know, this 
entails that E(g) is is an operator-valued analytic function in the sector 

Sk,s :={geC:0<\g\< rj{S); -{2K - 1)| + 5 < arg (g) < (2K + 3)| - 5} 

(ii) E(g,9) fulfills a strong asymptotic condition in 7Z (and thus, in particular, for 

oo 

g G Sk,s) and admits '^E i {Q 1 9)g l as asymptotic series; namely, there exist 

i=0 

A(5) > 0, C{5) > such that 



R N (g)\\ := \\E(g,9) - £ ^(0,%1 < AC^r((2X - 1)^/2)1^^ (2.10) 



E 

j=0 

as \g\ -> 0, (^,6>) eTZ, g e S KjS 



(iii) ^i(O,0) =E;(0,6), i = 0,l,..., #GR. 

Given the stability result (Assertion 2 of the present Theorem 1.3) the proof of (i) 
and (iii) is identical to that of [14], Lemma 2.5 (i) and is therefore omitted. We prove 
assertion (ii). Under the present conditions the Rayleigh-Schrodinger perturbation 
expansion is generated by inserting in (2.9) the (formal) expansion of the resolvent 
R Q (z,g,9):= [Q(g,0) - z]~ l : 

N-l 

R Q (z,g,9) = R Q (z,g,9) Y,[igWR v (z,0,0)] p + R Q (z)[igWR P (z,0,9)] N (2.11) 

p=0 

and performing the contour integration. Moreover (see once more [22], Section 5.7), to 
prove (2.10) it is enough to prove the analogous bound on D(g,9) and N(g,9). Since 
D(g) = P{0,9)P Q (g,9)P(0,t), we have, inserting (2.11) 

N-l 

D N (g,9) := D(g,9)- ]T A(0,% 1 

i=0 



P(0, 0)^- J r Rq(z, g, 9)[W(x)R v (z, 0, 9)] N P(0, , 
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By the analyticity and uniform boundedness of the resolvent Rq(z,q,9) in TZ (and 
hence in particular for g G Sk,s), it is enough to prove the estimate 

sup \\[igWR P (z,0,6)} N P \\ < AC N V{{2K - l)N/2)\g\ N (2.12) 

In turn, since P(O,0) = P(O,0), by the Combes-Thomas argument (see [22], Sect. 5 
for details) to prove (2.12) it is enough to to find a function / : R d — > R such that 

116^(0,0)11 < +oo; sup \W(x)e' f/N \ < (2.13) 

Now a basis in Ran(Pj) is given by m functions of the type 

Q(e^ Xl ,...,e^x d )e-^ 2 

where Q is a polynomial of degree at most m. Therefore both estimates are fulfilled by 
choosing / = a\x\ 2 with a = a(9) < 1/2. This condition is always satsfied if (g, 9) e TZ 
because |Im^| < it/ 4. This concludes the proof of the Theorem. 

Remark 

The summability statement just proved, called Borel summability for the sake of sim- 
plicity, is more precisely the Borel-Leroy summability of order q :— (K — l)/2. 

3 Conclusion 

Even though the object of main physical interest are the eigenvalues of H{g) rather than 
its singular values //&(<?) determined in this paper, the singular values yield a property 
that the eigenvalues cannot in general yield since the operator H(g) is not normal: 
namely, a diagonal form. If an operator is physically interesting a diagonalization of 
it is clearly useful. To examine this point in more detail, consider once again the 
canonical expansion (1.7) of Corollary 1.2: 

oo 

H(g)u = ]T fjL k (g) (u, ^ k )V^ k , u E D(H(g)) 

k=0 

Since both vector sequences {ipk} and {Vipk} are orthonormal we have 



{Vi/>k,H(g)iJ)i) = Hk(g)8k,i 



(3.14) 
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Moreover the orthonormal sequences {ipk} and {Vipk} are complete in the Hilbert 
space. Hence formula (3.14) is an actual diagonalization of H{g). The basis {ipk} 
acts in the domain, and the basis {Pipk} m the range. A complete diagonalization of 
the 'PT-symmetric but non-normal operator H(g) has been therefore obtained: the 
singular values fik(g) and the eigenvectors ipk (and thus also the vectors Vipk) are 
indeed uniquely defined by perturbation theory through the Borel summability. 
More precisely, the general formula (1.8) 

oo 

Hu = Vkiu, ipk)ip'k, u E D(H) 

fc=0 

which provides a diagonalization for an operator H with compact resolvent with respect 
to the pair of orthonormal bases {ipk} and {^}, requires a priori the computation of 
fi k and ipk as solutions of the spectral problem 

H*(g)H(g)^ = ^ (3.15) 

which represents an eigenvalue problem more complicated than H(g)<p = A0. The 
result of this paper means that the eigenvalue problem (3.15) can be replaced by the 
more tractable one 

H(g)ij> = iiVif> 

which can be solved by perturbation theory and Borel summability. 
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